On the dielectric susceptibility calculation in the incommensurate phase of K 2 Se04. 
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It is shown that the thermodynamic potential of the domain-like incommensurate (IC) phase of 
the K2Se04 crystal (viewed as a model for the IC-C transition) should be supplemented with a 
term, taking into account the local, Lorentz electric field. The latter qualitatively changes the result 
of calculation of the dielectric susceptibility for this IC structure by Nattermann and Trimper, J. 
Phys. C: Solid State Phys. 14, 1603, (1981), and gives phase transition to the ferroelectric IC phase 
obtained by Aslanyan, Phys. Rev. B 70, 024102, (2004). 



PACS numbers: 64.70.Rh, 64.60.-i 

In the present paper the behavior of the incommen- 
surate (IC) phase of the K^SeC^ crystal near the lock- 
in transition is discussed. We remind that K2Se04 un- 
dergoes on cooling two successive phase transitions: the 
IC transition at T$ = 129 K and the lock- in transition 
transition to the triple-period commensurate ferroelec- 
tric phase at T c = 93Km£ Nattermann and Trimper 3 
calculated the dielectric susceptibility \ of the crystal in 
the IC phase, and showed that \ diverges on approach- 
ing the lock- in transition point (i.e., a loss of stability 
of the IC phase with respect to the commensurate phase 
takes place). The dielectric susceptibility for the same 
IC structure was calculated also by Aslanyan^ and it 
was obtained that \ diverges not in the lock-in transi- 
tion point, but at a higher temperature, giving rise to the 
phase transition from the IC to the ferroelectric IC struc- 
ture. Spatial distributions of the IC modulation phase 
and polarization vector in the ferroelectric (or z- polar- 
ized) IC phased are depicted in Fig.l. Below we explain 
the origin of the difference between the results of the two 
calculations ^> and show that the thermodynamic poten- 
tial used for calculations by Nattermann and Trimper 3 
should be supplemented with a term, taking into account 
the local (Lorentz) microscopic electric field, induced by 
the dipoles (the neighboring pairs of the domains in Fig.l 
can be viewed as periodically arranged dipoles). In such 
a case, calculation of the dielectric susceptibility by Nat- 
termann and Trimper 3 should be revised, and it reduces 
to that carried out by Aslanyan^ 

The starting point in both calculations^ is the ther- 
modynamic potential of the following form: 
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In this expression (and below as well) the notations of 
paper [4] are used, where 770 and tp are, respectively, the 
amplitude and the phase of the IC modulation, k is the 
IC vector at the temperature of the IC transition (from 
the normal to the IC phase), P is the polarization vector 
in the z direction, and the IC modulation is given by 
tjq cos(kx + tp). 

However, potential (1) should be supplemented with a 



u(x) 



2ji 

4n 
3 

"T 



a) 



c) 



> — ' 






1 1 


. 25 I 


50 


1 1 



b) P 
r — ^ r — -\ r 



FIG. 1: a) Coordinate x dependence of the z -polarized IC 
structure's modulation phase u(x) = ip(x) + kx. Each plateau 
corresponds to the domain with + or - polarization vector P, 
which is along the ±z direction (curve b). Since the sizes of 
the adjacent domains are different, the integral polarization of 
any two adjacent domains (and, consequently, of the crystal 
as a whole) is not zero. Polarization of the crystal appears 
due to the displacement of the odd domain walls with respect 
to the even walls, which makes nonequal the sizes of adjacent 
domains. The coordinate x is given in multiples of the lattice 
parameter, c) Periodic disposition of the dipoles in the x — 
z plane. These dipoles correspond to the pairs of adjacent 
domains shown by the curve b). In the y and z directions 
the periods are equal to the crystal lattice parameters, while 
in the x direction the distance between dipoles equals to the 
width of the pair of domains. 



term of the following origin. The domain-like IC struc- 
ture, which develops from the sinusoidal one on cooling, 
can be polarized in the z direction via displacements of 
the odd domain walls with respect to the even walls in 
the x direction, as shown in Fig.l. The pairs of adjacent 
domains (which become non-equisize after such displace- 
ments) can be viewed as a periodically arranged dipoles 
as in Fig.l, c). However, any periodically (or chaoti- 
cally, as in the case of a gas) distributed dipoles, with 
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nonzero integral polarization, induce a microscopic elec- 
tric field E micr (the so called local Lorentz field) biassing 
each dipole. For the case of a gas or a cubic crystal, this 
field is equal to E m i cr — 4irP/3, where P is the spatially 
averaged polarization. The effect of the microscopic elec- 
tric field is implicitly taken into account for the suscep- 
tibility xo of the high temperature phase. However, in 
the IC phase, when the crystal acquires a new sublattice 
in form of the polarizing pairs of domains, one should 
additionally take it explicitly into account. So, for the 
orthorhombic domain lattice of K 2 Se04 the microscopic 
field is 



aP, withP 



where a is a coefficient of the same order as 47r/3 (in case 
of a cubic symmetry a = Att/3), and L is the crystal's 
length in the x direction. Due to the appearance of the 
field, biassing each dipole, the thermodynamic potential 
(f ) should be contributed by the additional term 
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which should be added to the expression under integral 
(1). ft should be emphasized a non-local character of this 
energy term (square of the integral) , which significantly 
changes the expected results of calculations. 

In what follows we show that the thermodynamic po- 
tential of the crystal with the domain-induced micro- 
scopic electric field can be reduced to the potential, ana- 
lyzed in the paper [4] , and the behavior of the system in 
such a case should be that described ibidem. Minimizing 
the potential over P{R) and P, under the condition 



(P(r) - P)dx = 0, 



which should be added to the potential equation with the 
Lagrange multiplier, one obtains 
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where u = kx + <p(x). One can substitute P(R) and P 
back into the potential expression and obtain: 



16 



4L 2 ( 

We once more note a non-local character of the coupling 
in this potential, which is given by the square of the in- 
tegral in the equation. Such a coupling is induced by the 



above-introduced microscopic electric field, and it disap- 
pears at a = 0. Eq.(2) allows one (via minimization over 
f{R)) to derive the equation for the temperature evolving 
of the IC structure. For the case a = the solution of the 
corresponding equation is well studied^ and it gives IC 
domains of equal sizes, which are increasing (logarithmi- 
cally diverging) on cooling towards the lock-in transition. 
The dielectric susceptibility for such a behavior was cal- 
culated by Nattermann and Trimper^ who showed that 
it diverges at the lock-in transition point. However for 
the case a ^ 0, eq.(2) significantly differs from that for 
a = 0. It is easy to check that eq.(2) is equivalent to: 

3 = J dx{^ V 2 (^-) 2 + ^r 1 e cos(6kx + 6v) + 



Piri cos(3kx + 3<p) + ^-Pj 2 }, (3) 



where Pi is a spatially constant vector and 
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In order to check that at a ^ eq.(3) should give the 
same solutions for the phase ip(R) as eq.(2), it is sufficient 
to minimize eq.(3) over P 1; and, as a result, return to 
cq.(2). 

On the other hand, eq.(3) coincides by form with 
eq.(l), with P replaced by the spatially constant vector 
Pi. Just the potential of eq.(3) with constant Pi was used 
by Aslanyan4 (though the approach was not commented 
and sufficiently clarified there) for checking the stability 
of the IC phase, and it was shown that within the IC 
structure a phase transition should take place from the 
IC to the ferroelectric IC phase, prior to the lock-in tran- 
sition on cooling. The spatial distribution of the phase 
u = ip + kx and polarization P(P) in this ferroelectric 
phase are those shown in Fig.li4 

In other words, for calculation of the dielectric suscep- 
tibility x 1 = d 2 &/dP 2 at P = one should minimize 
potential (1) with respect to (p, and find the solution of 
the equation: 
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In the paper ref. [4] the corresponding solution was found 
for the case P = const in this equation (which corre- 
sponds to minimization of eq. (3)), and in ref. [3] that 
was obtained for a spatially dependent P(R). 

So, we showed that thermodynamic potential (1), used 
for calculation of the dielectric susceptibility, should be 
supplemented by the term, contributed from the domain- 
induced Lorentz local electric field. As a result, one can 
see that calculations of the dielectric susceptibility by 
Nattermann and Trimper— give correct result only for a 
pair of the adjacent domains, which is separated from 
the crystal (i.e. for a single dipole). For a large number 
of pairs of domains, periodically disposed in space as in 
Fig.l (in the real IC phase), the coupling between the 
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dipoles should be taken into account, which give phase 
transition to the ferroelectric IC phase with non-equisizc 
IC domains. This coupling gives contribution of a neg- 
ative sign to the potential, and makes ferroelectric IC 
phase more favorable compared to the equisize domain 
IC structure on approaching the lock-in transition. 

As an experimental confirmation of the presented the- 
ory may be pointed out the second-harmonic generation 
in the IC phases of the K 2 SeO/, (NH 4 ) 2 BeF^ and, un- 
der some assumptions, of the quartz^ crystals, which 
is discussed in ref.[4], and the present paper should be 



viewed as an additional foundation of its evidence. It is 
also worth to note that the developed theory allows one 
to explain such a surprising observations as a significant 
drop in the K2Se04 elastic constant C55 (which corre- 
sponds to the transversal u z sound propagating in the x 
direction) near the lock-in transition point^, and obser- 
vation of a depolarized, overdamped Raman scattering 
in the z(yz)x geometry^ in the IC phase. These obser- 
vations unambiguously follow from the presented theory, 
and will be published elsewhere. 
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